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Abstract. We used our database of VUT-UVES quasar spectra to build up a sample of 70 Damped Lyman-a (DLA) or strong 
sub-DLA systems with total neutral hydrogen column densities of log ^(111) i 20 and redshifts in the range 1.7 < z-^bf. < 4.3. 
For each of the systems, we measured in an homogeneous manner the metallicities relative to Solar, [X/H] (with X = Zn, or 
S or Si), and the velocity widths of low-ionization line profiles, AV. We provide for the first time evidence for a correlation 
between DLA metallicity and line profile velocity width , which is detected at the 6 .10" significance level. This confirms the 
trend previously observed in a much smaller sample bv IWolfe & Prochaskal Jl998h . The best-fit linear relation is [X/H] = 
1.55(±0. 12) log AV - 4.33(±0.23) with AV expressed in km s"'. The slope of the DLA velocity-metallicity relation is the same 
within uncertainties between the higher (z^bs > 2.43) and the lower (z;,bs ^ 2.43) redshift halves of our sample. However, the two 
populations of systems are statistically different. There is a strong redshift evolution in the sense that the mean metallicity and 
mean velocity width increase with decreasing redshift. We argue that the existence of a DLA velocity-metallicity correlation, 
over more than a factor of 100 spread in metallicity, is probably the consequence of an underlying mass-metallicity relation 
for the galaxies responsible for DLA absorption lines. Assuming a simple linear scaling of the galaxy luminosity with the 
mass of the dark-matter halo, we find that the slope of the DLA velocity-metallicity relation is consistent with that of the 
luminosity-metallicity relation derived for local galaxies. If the galaxy dynamical mass is indeed the dominant factor setting 
up the observed DLA velocity-metallicity correlation, then the DLA systems exhibiting the lowest metallicities among the 
DLA population should, on average, be associated with galaxies of lower masses (e.g., gas-rich dwarf galaxies). In turn, these 
galaxies should have the lowest luminosities among the DLA galaxy population. This could explain the difficulties of detecting 
high-redshift DLA galaxies in emission. 
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1. Introduction 

Over the past decade, significant progress in our understanding 
of early galaxy evolution has been made with large samples 
of high-redshift galaxies drawn from deep multi-band imag- 
ing Isteidel et al. 2003, and references therein). However, even 

* Based on data gathered at the European Southern Observatory 
(ESO) using the Ultraviolet and Visible Echelle Spectrograph (UVES) 
installed at the Very Large Telescope (VLT), Unit-2, Kueyen, on Cerro 
Paranal in Chile. 



before the first surveys for Lyman-Break Galaxies (LBGs) 
had begun, samples of DLA absorbers observed on the lines- 
of-sight to d istant quasars had been constructed llWolfe et alJ 
Il986l I1995I) . These absorbers were thought at the time to be 
the best carriers of information on the population of high- 
redshift galaxies, but, despite many attempts to identify the 
galaxies responsible for DLA absorption lines [hereafter called 
DLA galaxies], on ly very fe w could be detected i n emis- 
sion (see, e.g.. M0 ller & Warren 1993: Diorg ovski et al. 1996: 
iLeibundgut & Robertsonll999HFvnbo et al.ll999l:l]V[0ller et alJ 
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l2002i l2Q04 IChristensen etal]l20Q4l: IWeatherlev et^hOQ^ . 
However, there is little doubt that DLA systems arise from 
the densest regions of the Universe and are closely associated 
with galaxies. It is therefore crucial to establish the connection 
between the absorption-selected DLA systems and emission- 
selected galaxies. In addition, the detailed information that be- 
comes available only through the combination of morphology, 
colour and luminosity, with QSO absorption-line spectroscopy, 
makes these galaxy/absorber associations unique laboratories 
to study t he phys ical processes at work during galaxy forma- 
tion (see lPei et alJl999l) . 

Progress in this field has been slow. Firstly, a huge amount 
of work is needed to derive important parameters in DLA sys- 
tems such as gas kine matics, metallicity, or dust and molecu- 
lar fra ctions (see, e.g., Pettini et al. '"1997I 'Prochaska & Wolfe' 
19971 [l998. 1999; Prochaska et al. 2001; Ledoux etal. 1998, 



20031) . Secondly, as mentioned above, the known high-redshift 



DLA systems have proved to be very difficult to detect in emis- 
sion. This has caused some confusion and for a while sug- 
gestions were put forward that DLA ab sorbers may no t be 
relat ed to high-redshift galaxies at all. Mo et alJ ( Il998l ) and 
iHaeh nelt et al. ( 1998b, 2000) resolved this issue showing that 
the difficulty of detecting high-redshift DLA systems in emis- 
sion is an unavoidable consequence of the absorption cross- 
section selection which tends to reveal faint galaxies because 
they have an inte grated cross-sectio n larger than that of bright 
galaxies (see also Fvnbo et al.* 1 9991) . 

Recently, M0ller et al. ( 2004,) tentatively suggested that, if 
a galaxy luminosity-m etallicity relation similar to that observed 
at ^ z ^ 1 (e.g. lGarn ett"2002^. 'Kobulni ckv et"ai] l2003t 
iLamareille et alJ [20041 iTre monti et al. 2004) was already in 
place at high redshifts, then it would be possible to significantly 
increase the DLA galaxy detection probability by carefully se- 
lecting DLA systems with the highest metallicities. In fact, the 
few DLA galaxies that have to date been identified in emission 
do give support to the conjecture that a luminosity-metallicity 
relation was already in place at z x 2-3, altho ugh the re- 
sult i s only marginally statistically significant ( M0ller et alJ 
l2004 see also Christensen et al., in prep.). This is in line 
with the near-Solar or even super-Solar metallicities derived 
for bright Lyman-break or bright K- band selected galaxie s at 
similar redshifts (.Shaolev et aL . 2004ilde Mello et alJl2004 . A 
mass-metallicity relation has recently been put into evidence 
for UV -selected star-forming galaxies at z ~ 2.3 bv lErb et alJ 

toodi . 

In this paper, we provide for the first time evidence for the 
existence of a velocity-metallicity correlation for high-redshift 
DLA galaxies that could be the consequence of an underly- 
ing mass-metallicity relation for the galaxies responsible for 
DLA absorption lines. From the observation of a sample of 
17 DLA systems at z^bs < 3. IWolfe & Prochaskal ( fl998 ) pre- 
viously showed that the DLA systems exhibiting the largest 
line profile velocity widths span a narrow range of high metal- 
licities. However, these authors also suggested that systems 
with sma ll velocity widths sp an a wide range of metallicities. 
Recentlv. lPeroux et alJ fe003) found a hint of an increase of the 
mean DLA metallicity with increasing velocity width, but the 
statistical significance of their result is low. In this paper, we 



use our database of VLT-UVES quasar spectra to build up a 
sample of 70 DLA or strong sub-DLA systems with total neu- 
tral hydrogen column densities of log A^(H i) ^ 20 and redshifts 
in the range 1.7 < Zabs < 4.3. We present new, homogeneous 
measurements of DLA metallicities and line profile velocity 
widths in Sect.|5]and the observed velocity-metallicity relation 
in Sect. 13 We discuss in Sect.|3]the use of the DLA gas kine- 
matics as a proxy for the mass of DLA galaxies and the possi- 
bility of the existence of a mass-metallicity relation for high-z 
DLA galaxies. We conclude in Sect.lS] 

2. UVES DLA sample 

2.1. Metallicity measurements 

Most of the systems in our sample were select ed from the 
follow-up of the Large Bright QSO Survey (LBOS: IWolfe et alJ 
,1995.) and observed at the ESO VLT with UVES between 2000 
and 2004 in the cour se of a systematic search for molecular hy - 
drogenatzabs > 1-8 (Pet itiean et al . 2000; Ledoux et al."20Ii2|). 
Our total sample comprises 57 DLA systems (logA^(Hi) > 
20.3) and 13 strong sub-DLA systems with total neutral hy- 
drogen column densities in the range 20.0 ^ log A^(H i) < 20.3. 
This is only slightly lower than the classical definition for DLA 
systems to ensure that these absorbers are mostly neutral and 
share the sam e physical nature as classical DLA systems (see 
IViegaJl995l). 

We have carefully measured or remeasured total neutral hy- 
drogen column densities, log A^(H i), and average DLA metal- 
licities for all the systems in our sample. Results are summa- 
rized in Table [J The absorption line analysis was performed 
in an homogeneous manner using standard Voigt-profile fit- 
ting techniques adopting the oscillator strengths compiled by 
Morton (2003). For the damping coeffi cients, we also adopted 
here the compilation by MortonI (|2P03), which results in some 
cases in a sl ight increase of logA^(Hi) values compared to 
iLedoux et alJ J2003). Total metal column densities were de- 
rived as the sum of the column densities measured in in- 
dividual components of the line profiles. Average gaseous 
metallicities relative to Solar, [X/H] = log[A^(X)/A?(H)] - 
log[A^(X)/A^(H)1p, were calculated using Solar abundances 
listed in Morto'J(l2003l) which are based on meteoritic data 
from iGrevesse & Sauvall 1I2OO2I) . To avoid problems related 
to possible depletion onto dust grains, metallicities given in 
Tabled were computed for elements that are known to deplete 
very little in the ISM of the Galaxy. The reference element was 
taken to be X ^Zn when Z n 11 is detected, or else either S or Si 
was used (see lLedoux et al..2003.. for a discussion). 

A noticeable property of this large dataset is that it sam- 
ples well both ends of the DLA metallicity distribution, from 
[X/H] ^ -2.6 up to about half of Solar (see TableHJ. 



2.2. Low-ionization line kinematics 

For each of the systems, we also determined the velocity 
widths of metal absorption line profiles. Low-ionization tran- 
sition lines that are not strongly saturated were selected to 
trace the kinematics induced predominantly by gravity. For 
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Fig. 1. Measurements of the velocity width of low-ionization line profiles, AV. The different panels show the apparent optical 
depth in the selected transition line for a given system: Siii/11526 at Zabs = 2.418 toward Q 0112-306 (upper left), Feii/11608 
at zabs = 3.025 toward Q 0347-383 (upper right). Si ii^l526 at Zabs = 2.622 toward Q 0405-443 (lower left), and Fe 1^1608 at 
Zab,s - 2.426 toward Q 2348-011 (lower right). The top part of each panel shows the integration of the apparent optical depth 
starting at wavelength Ai. The horizontal bar below that corresponds to the measurement of AV when 5% of the total apparent 
optical depth is avoided at both edges of the profile. 



high-ionization lines, the velocity widths could be dominated 
by peculiar ejection of hot gas. We measured the line profile 
velocity widths follow ing [Prochaska & Wolfe ( 1997; see also 
lHaehneltetal]ll998bh . We used the criterion that the residual 
intensity / of the strongest absorption feature in the selected 
line profile must satisfy 0.1 < I/Ic < 0.6, where is the in- 
tensity level of the adjacent continuum. This criterion selects 
transitions that are neither strongly saturated (in which case the 
optical depth cannot be derived and the velocity width could be 
overestimated), nor too weak (in which case the velocity width 
could be underestimated because part of the gas would be un- 
detected). For the few systems for which none of the observed 
lines satisfies the above criterion, we used the mean value of 
the velocity widths calculated from (i) a line slightly more sat- 
urated, and (ii) a line slightly weaker, than what the criterion 
specifies (see TableQJ. From visual inspection of the strongest 



low-ionization line profiles of a given system, we established 
the velocity range over which the previously selected line pro- 
file should be integrated (this corresponds to the interval [Ai ;A2] 
in Fig. [3. 



In order to ease com parison with previous works (e.g., 
IProchaska & W olfe 1997), we then calculated the line profile 
velocity width, AV, as c[A(95%)-A(5%)]/Aa, where A(5%) and 
A{95%) are the wavelengths corresponding to, respectively, the 
five per cent and 95 per cent percentiles of the apparent optical 
depth distribution, and Aq is the first moment of this distribu- 
tion (see Fig. [3. Note that all our spectra have a typical signal- 
to-noise ratio larger than 30. Excluding the extended wings of 
the line profiles avoids taking into account satellite components 
that are not strictly related to the bulk of the systems. 
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Fig. 2. Average metallicity of each of the 70 DLA or strong 
sub-DLA systems in our sample, [X/H], vs. the velocity width 
of their low-ionization line profiles, AV, displayed on a loga- 
rithmic scale. A positive correlation between the two quanti- 
ties is detected at the 6.1cr significance level using a Kendall 
rank correlation test. The linear least-square bisector fit is 
shown as a solid line. Sub-DLA systems, with total neutral 
hydrogen column densities 20 5 logA^(Hi) < 20.3, are indi- 
cated by filled triangles. There are a few DLA systems with 
c|Zem -2absl/(l +Zabs) ^ 5000 km s"' . For the sake of complete- 
ness, they are shown here as empty circles but they were not 
considered in the analysis. 

3. The velocity-metallicity relation 

3.1. Detected correlation 

In Fig. 12 we plot on a logarithmic scale the average metallic- 
ity of each of the 70 DLA or strong sub-DLA systems in our 
sample, [X/H], versus the velocity width of their low-ionization 
line profiles, AV. A positive correlation between the two quan- 
tities is detected at the 6Acr significance level using a Kendall 
rank correlation test. Note that the DLA and strong sub-DLA 
populations are statistically indistinguishable even though the 
mean metallicity is slightly larger among sub-DLA systems. 
The Kolmogorov-Smirnov test probability that the two ve- 
locity width distributions (resp. metallicity distributions) are 
drawn from the same parent population is Pks = 0.91 (resp. 
^'ks - 0.72) in the two-sided case. In the following, we there- 
fore consider the DLA and strong sub-DLA systems in our 
sample altogether 

The typical measurement uncertainties in velocity width 
(+0.02 dex) and metallicity (+0.10 dex) are small compared 
to the intrinsic scatter of the data points (see Fig.|2}. We thus 
fitted the data us ing the linear least-square bisector method 
lllsobe et al.l ll990). For the correlation fits, we did not include 
the DLA systems with c|zem - ■Zabsl/(1 + Zabs) ^ 5000 km 
s where Zem is the QSO emission redshift. However, these 
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Fig. 3. Histogram of the projected position, t/pioj, of each sys- 
tem along the best-fit correlation relation plotted as a solid line 
in Fig. 12 with t/pi-oj - (resp. c/pi-oj = 1) at the intersection of 
this line with the horizontal line corresponding to [X/H] - -3 
(resp. [X/H] = 0). The shaded (resp. hashed) histogram repre- 
sents the higher (resp. lower) redshift half of the sample. 

Zabs ~ Zem systcms are not associated with the central engine 
of the quasar, nor ejected by the quasar, but rather associated 
with d ense gas in its surroundings (e.g..'Petitiean et al."l 994t 
ISrianand & Petitiean 2000) . It can be seen in Fig. |2] that their 
inclusion would not affect the observed correlations and that 
their metallicities and line profile velocity widths are represen- 
tative of those of the overall DLA population. 
The best-fit linear relation is: 

[X/H] = L55(+0.12)logAy-4.33(±0.23) (1) 

with Ay expressed in km s"' . 

The observed high-redshift velocity-metalUcit y correlation 
presented here is consistent with th e results of 'Ne stor et alJ 
( 2003.) (see also lTurnshek et al.l2005l) . Using the Sloan Digital 
Sky Survey (SDSS), these authors found that for strong low- 
ionization Mg II absorbers at 1 $ Zabs ~ 2 the Mg ii/12796 equiv- 
alent width, or equivalently the absorption line velocity spread, 
is correlated with the metallicity. In addition, they showed that, 
within the large equivalent width regime, the average metalUc- 
ity is larger at lower redshift. 

3.2. Redshift evolution 

The median redshift of our sample is Zmed = 2.43. It can be 
shown that the two sub-samples of systems with, respectively, 
Zabs > 2.43 and Zabs ^ 2.43 differ significantly. The median 
DLA metallicity and median DLA velocity width increase with 
decreasing redshift: [X/H] - -1.59 and AV - 69 km s"', re- 
spectively, [X/H] - -1.15 and AV = 92 km s in the higher, 
respectively, the lower redshift half of the sample. This resem- 
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Fig. 4. Same as Fig.|2lbut for the sub-sample of systems with 1.7 < Zabs < 2.43 {left hand-side panel) or the sub-sample of 
systems with 2.43 < Zabs < 4.3 {right hand-side panel). Positive correlations between the two quantities are detected in both the 
lower and the higher redshift sub-samples with significance levels of 4.6cr and 3.9cr, respectively. Linear least-square bisector 
fits to each sub-sample are shown as solid lines in the corresponding panels; the dotted lines are the fits to the other sub-sample. 
There are a few DLA systems with c|zem - Zabsl/(1 + ^abs) ^ 5000 km s"' in our sample. For the sake of completeness, they are 
shown here as empty circles but they were not considered in the analysis. 



bles the point bv lWolfe & Prochaskal ll 19981) that the kinemat- 
ics and metallicities of the Zabs > 3 and Zabs < 3 DLA samples 
could show significant differences. 

To investigate this behaviour further, we have calculated 
for each system in our sample the projected position, t/pioj, 
along the best-fit correlation relation ([X/H] vs. AV) derived 
in Subsect. [3.1l and drawn as a solid line in Fig.|2 We plot in 
Fig-E]the histograms of t/pioj for the two redshift sub-samples. It 
is apparent that the two histograms differ in the sense that t/pioj 
is larger at lower redshift. This is expected as dpiq increases 
with increasing velocity width and increasing metallicity. A 
Kolmogorov-Smirnov test confirms this: the two populations 
have only 1 % chance to be drawn from the same parent popu- 
lation. This difference is also apparent in Fig. 0] where we plot 
[X/H] versus AV for the two redshift sub-samples separately. 

Considering the two redshift sub-samples separately, the 
best-fit linear relations are: 

[X/H] = 1.45(+0.17)logAy-4.01(+0.33) (2) 
for 1.7 < Zabs < 2.43, and: 

[X/H] = 1.47(+0.17)logAy-4.32(+0.32) (3) 

for 2.43 < Zabs < 4.3, with AV expressed in km s"'. 

The significance levels of the correlations in the lower, re- 
spectively, the higher redshift half of the sample are A.6cr and 
3.9cr, respectively (see Fig.|3. In addition, the Pearson corre- 
lation coefficients are r - 0.72 and 0.63, respectively, showing 
that even in the high-redshift sub-sample a linear relation is a 
fairly good description of the data. It is very interesting to note 
that the correlation relations do not change significantly with 



redshift although there is a statistically significant increase in 
both [X/H] and AV with decreasing redshift. This is discussed 
in Sect.m 

3.3. Observational scatter 

Scatter in the data points is expected due to random impact pa- 
rameters through the absorbing galaxies and, indeed, the scat- 
ter observed in Fig.|3]is much larger than the metallicity mea- 
surement uncertainties. Negative radial gradients in metallicity 
like those observed in the discs of large nearby spirals could 
also contribute to the scatter of the data points. T his effect is 
probably not very important however. For instance. Ichen et all 
( 2005) derived from a sample of six z < 0.65 galaxy-DLA pairs 
a metallicity gradient of only -0.041 + 0.012 dex per kpc from 
galactic centre to 30 /z ' kpc radius. In addition, the magnitude 
of such gradient s in the discs of nea rby spirals has been ques- 
tioned recently jSresolin et alJ2004l) . 

We note that there are in our sample a few systems de- 
parting from the general trend. The system at Zabs = 1.973 
toward Q 0013-004 has low-ionization hnes extending ove r 
up to ~1100 km s"' in velocity space (Petitjea n et aljEoo3) . 
which is much larger than for the other systems in our sam- 
ple. The DLA system at Zabs = 2.622 towai'd Q 0405-443 has 
a low metallicity, [Si/H] = -2.04, and a comparatively large 
velocity width for this metallicity, AV = 182 km s"' (see 
Fig H lower left-hand side panel). An inspection of these two 
cases indicates that the low-ionization line profiles are clumpy, 
being made of four and two well-separated clumps, respec- 
tively. These systems could arise in galaxy groups, tidal stream- 
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ing in galaxy mergers or forming galactic structures (see, e.g., 4.2. Implications and prospects 
' "^^^ ilJ I 19981: ■ ' 



Ledoux et a 



Mailer et al 



"l998"; N ulsen etalJI 19981: iHaehnelt et ; 
12001) . The third peculiar system, the sub-DLA 



^ Il998bt 



system at Zabs - 2.287 toward Q 2332-094, is the highest 
metallicity absorber in our sample with [Zn/H] = -0.33. It has 
a surprisingly small velocity width for such a high metallicity 
however: its profile is a blend of two sharp metal lines resulting 
in a velocity width of AV = 43 km s"' . It is not unexpected to 
find such deviant cases in a sample that large. 



4. Discussion 

4.1. Kinematics as a proxy for the mass of DLA 
galaxies 

The existence of a DLA velocity-metallicity correlation, over 
more than a factor of 100 spread in metallicity, can be under- 
stood as the consequence of an underlying mass-metallicity re- 
lation for the galaxies responsible for DLA absorption lines. 
Peculiar ejection of hot gas should indeed primarily affect the 
kinematics of high-ionization lines such as C iv and Si iv. Our 
measurements of profile velocity widths in DLA systems are 
based on low-ionization lines which should instead be dom- 
inated by motions on galactic scale governed or induced by 
gravity. 

For disc galaxies, the rotation velocity is a direct mea- 
sure of the galaxy's total mass. Lines-of-sight that do not pass 
through the centre of the galaxy will not trace the full depth 
of the potential well and, therefore, will tend to show smaller 
veloc ity dispersions than the ro tat ion velocity (see the mod - 
els bv lProchaska & Wolfell997l and lWolfe & Procha ska 1998). 
Hence, for random lines-of-sight through a large sample of disc 
galaxies, there should be a mean relationship between mass and 
profile velocity width, albeit with a large scatter induced by the 
range of impact parameters and inclination angles probed by 
the observations. 

Infall/outflow of gas, or merging of galaxy sub-clumps, will 
also produce kinematically broadened line profiles with veloc- 
ity widths scaling as the infall/outflow velocities, which again 
scale as the depth of the combined potential well of the galax- 
ies or mergers. Simulations have shown that in that case the 
line profile velocity width is a good indicator of the circular 
velocity of the underlying dark-matter halo with AV ~ 0.6v p 
jNulsen etalJ 'l998: H aehnelt et a l. 1998b: Ma iler et al.l2 001). 
The scatter in this relation is about a factor of two and corre- 
sponds approximately to the width of the correlation shown in 
Fig.E 

In addition, there is a positive correlation between the pro- 
jected stellar mass density and the neutral hydrogen column 
density of DLA systems and a good correspondence in the spa- 
tial distribution of stars and DLA systems in simulations in- 
cluding star forma tion, supernova fee dback and feedback by 
galactic winds (see lNagamine et al'l2004.) . 



For reasons discussed above, we assume here that the dynami- 
cal mass of DLA galaxies is the dominant factor setting up the 
observed DLA velocity-metallicity correlation. 

4.2.1. Mass-metallicity relation 

The slope of the DLA velocity-metallicity relation is the same 
within uncertainties for the two redshift halves of our sample 
(see Subsect. [3.2> . There is a possible increase of the intercept 
point of this relation with decreasing redshift (see Fig. |3 but 
this result is not statistically significant due to the large scatter 
of the data points around the mean relations (see Eqs.|2land|3}. 
However, the two redshift sub-samples differ in the sense that 
the median DLA metallicity and median DLA velocity width 
increase with decreasing redshift (see Subsect. f3.2l and Fig.|3. 
This suggests that galaxy halos of a given mass (resp. a given 
metallicity) are becoming more metal-rich (resp. less massive) 
with time. 

This is consistent with the results o f lSavaglioetaTI(l200.4 
who proposed a redshift-dependent galaxy mass-metallicity re- 
lation from the study of 0.4 < z < 1.0 galaxies selected from 
the Gemini Deep Deep Survey and the Canada-France Redshift 
Survey. Note also that a mass-metallicity relation has recently 
been put into evidence at z ~ 2.3 for UV-selected star-forming 
galaxies (.Erb et al„2006.) . 

4.2.2. Luminosity-metallicity relation 

From C old Dark Matter (CDM) simulations, iHaehnelt et alJ 
(Il998l j) have shown that the velocity width of DLA systems, 
AV, can be related statistically to the circular velocity of the 
underlying dark-matter halo, = (GM/ry^^, where M is the 
mass in a sphere overdense by a factor of 200 compared to the 
mean cosmic density. They found AV ~ 0.6v'c. According to 
[Haehnelt et alJ (ll998all200Cll) . the luminosity function of z ~ 3 
galaxies can be reproduced if a simple linear scaling of the lu- 
minosity with the mass of the dark-matter halo is assumed, i.e., 
Mr = -7.5 log(vc/200 km s"') + m^, where mr is the galaxy 
apparent /?-band magnitude and = 26.6 for the A-CDM 
model. Using our best-fit to the velocity-metallicity relation for 
1-7 < Zabs ^ 2.43 DLA systems (Eq.|2jl, and Zmed = 2.09 for 
this sub-sample, we derive: 

[X/H] = -0.19(+0.02)MR-4.76(±0.42) + 0.19(±0.02)/iOj (4) 

where Mr is the galaxy absolute /?-band magnitude and Kr 
is the /T-correction in the 7?-band. Using our best-fit to the 
velocity-metallicity relation for 2.43 < Zabs < 4.3 DLA sys- 
tems (Eq.|3} leads to a similar result. 

It is striking to note that the slope of this DLA 
lu minosity-metalli c ity re lation is consistent with that derived 
by iTremonti et alJ (|2004) for the luminosity-metallicity rela- 
tion for z ~ 0.1 galaxies selected from the SDSS, [O/H] = 
-0.185(+0.001)Mb - 3.452(+0.018). The correction from the 
/?-band at high redshift to the Z?-band at low redshift cor- 
responds to an additi onal factor that co rresponds to a non- 
positive /T-correction jKinnev etalJl996l) in Eq. |3 Therefore, 
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the intercept points of the two luminosity-metallicity relations 
are different in the sense that galaxies of a given luminosity 
(resp. a given metallicity) are becoming more metal-rich (resp. 
fainter) with time. 

5. Conclusions 

Using a sample of 70 DLA or strong sub-DLA systems with 
total neutral hydrogen column densities of log A^(H i) ^ 20 and 
redshifts in the range 1 .7 < z^bs < 4.3, we have shown that there 
is a correlation between metallicity ([X/H]) and line profile ve- 
locity width (AV) at the 6.1cr significance level. The best-fit 
linear relation is [X/H] = 1.55(+0.12) log AV - 4.33(+0.23) 
with Ay expressed in km s~'. We argued that the existence of 
a DLA velocity-metallicity correlation, over more than a factor 
of 100 spread in metallicity, is probably the consequence of an 
underlying mass-metallicity relation for the galaxies responsi- 
ble for DLA absorption lines. Assuming a simple linear scaling 
of the galaxy luminosity with the mass of the dark-matter halo, 
we found that the slope of the DLA velocity-metallicity rela- 
tion is consistent with that of the luminosity-metallicity rela- 
tion derived for local galaxies. If the galaxy dynamical mass 
is indeed the dominant factor setting up the observed DLA 
velocity-metallicity correlation, then the DLA systems exhibit- 
ing the lowest metallicities among the DLA population should, 
on average, be associated with galaxies of lower masses. 

Eq. 0] implies that the more than two orders of magni- 
tude spread in DLA metallicity could reflect a more than ten 
magnitudes spread in DLA galaxy luminosity. Even though 
low-mass galaxies, i.e., gas-rich dwarf galaxies, can undergo 
periods of intense star formation activity and, in this case, 
have high luminosities in the UV, it is a fact that, on aver- 
age, they show lo wer star-format ion rates than more massive 
galaxies ( Brinchm ann et al J2004t see also Okoshi et al. 2004). 
This may well explain the diflicu lty of detectin g high-redshift 
DLA galaxies in emission (e.g., iKulkarni et al...200a .2001) . 
Furthermore, the non-detection of Lya emission from the 
galaxies responsible for low-metallicity DLA systems, down 
to Lya fluxes fainter than most of the Lyo' emitters from the 
deep survey of Fvnbo et al. (2003), could be a consequence of 
their low masses and their correspondingly, on average, low 
star-formation activity. 

A significant fraction of DLA galaxies could be actively 
forming stars for some period of time, but due to their small 
masses the bursts of star-formation would not be powerful 
enough and/or would have too short life-times. Lya emission 
would simply be too faint or would be strong for a too short 
spell of time to be detected by the current generation of 8- 10 m 
class telescopes. Conversely, the existence of a DLA mass- 
metallicity re lation can explain the recent, tentative result by 
iMoUer et all (12004 ) that the few DLA systems with detected 
Lya emission have higher than average metallicities (see also 
Christensen et al., in prep.). This should be confirmed by ad- 
ditional deep imaging of the fields of QSOs with carefully se- 
lected DLA absorbers. 

The DLA velocity-metallicity correlation relation studied 
in this paper also needs to be investigated in the context of new 
high-resolution smoothed-particle hydrodynamics simulations 



including the eff'ects of feedbac k in a self-consistent manner 
(see, e.g.. lNagamine et alJ2004l) . 
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Table 1. UVES DLA sample: average metallicities and velocity widths of low-ionisation line profiles * 



Quasar 


Other name 




2abs 


logiV(Hi) 


[X/U]' 


X 


Ay 

(km s-i) 


Selected 
transition lines ^ 


Q 0000-263 


LBQS 0000-2619 


4.11 


3.390 


21.40 + 0.08 


-2.06 + 0.09 


Zn 


33 


Oa950 


Q 0010-002 


LBQS 0010-0012 


2.15 


2.025 


20.95 + 0.10 


-1.43 + 0.11 


Zn 


32 


Siiai808 


Q 0013-004 


LBQS 0013-0029 


2.09 


1.973 


20.83 ± 0.05 


-0.59 ± 0.05 


Zn 


720 


Fe 11^1608 


Q 0058-292 


LBQS 0058-2914 


3.09 


2.671 


21.10 + 0.10 


-1.53+0.10 


Zn 


34 


Si nil 808 


QOlOO+130 


LBQS 0100+1300 


2.69 


2.309 


21.35 + 0.08 


-1.58 + 0.08 


Zn 


37 


CriU2056 


Q 0102-190 


LBQS 0102-1902 


3.04 


2.370 


21.00 + 0.08 


-1.90 + 0.08 


S 


17 


Fe 11^2374 


Q0102-190 


LBQS 0102- 1902 


3.04 


2.926 


20.00 + 0.10 


-1.50 + 0.10 


Si 


146 


Si 11/11526 


Q 01 12-306 


LBQS 01 12-3041 


2.99 


2.418 


20.50 ± 0.08 


-2.42 + 0.08 


Si 


31 


Siiai526 


Q 01 12-306 


LBQS 01 12-3041 


2.99 


2.702 


20.30 + 0.10 


-0.49 + 0.11 


Si 


218 


Fe nil 608 


Q 01 12+030 


LBQS 01 12+0300 


2.81 


2.423 


20.90 + 0.10 


-1.31+0.11 


S 


112 


Fe Hi 1608 


Q 0135-273 


CTS 1005 


3.21 


2.107 


20.30 + 0.15 


-1.12 + 0.16 


S 


82/103 


Feii/12260/Fen/12586 


Q 0135-273 


CTS 1005 


3.21 


2.800 


21.00 + 0.10 


-1.40 + 0.10 


S 


65/67 


SiLil259/Fenil608 


Q 0216+080 


LBQS 0216+0803 


2.99 


1.769 


20.30 + 0.10 


-1.00 + 0.12 


Zn 


75 


Fe 11^2374 


Q 0216+080 


LBQS 0216+0803 


2.99 


2.293 


20.50 + 0.10 


-0.70 + 0.11 


Zn 


104 


Siiai808 


Q 0336-017 


LBQS 0336-0142 


3.20 


3.062 


21.10 + 0.10 


-1.41+0.10 


Si 


67 


Si nil 808 


Q 0347-383 


LBQS 0347-3819 


3.22 


3.025 


20.73 ± 0.05 


-1.17 + 0.07 


Zn 


93 


Fe nil 608 


Q 0405-443 


CTS 247 


3.02 


1.913 


20.80 + 0.10 


-1.03 + 0.10 


Zn 


98 


Si nil 808 


Q 0405-443 


CTS 247 


3.02 


2.550 


21.15 + 0.15 


-1.36 + 0.16 


Zn 


165 


Si nil 808 


Q 0405-443 


CTS 247 


3.02 


2.595 


21.05 + 0.10 


-1.12 + 0.10 


Zn 


79 


Si nil 808 


Q 0405-443 


CTS 247 


3.02 


2.622 


20.45 + 0.10 


-2.04 + 0.10 


Si 


182 


Sinil526 


Q 0450-131 




2.25 


2.067 


20.50 ± 0.07 


-1.62 + 0.08 


S 


148 


Fenil608 


Q 0458-020 


LBQS 0458-0203 


2.29 


2.040 


21.70 + 0.10 


-1.22 + 0.10 


Zn 


88 


Crni2056 


Q 0528-250 


LBQS 0528-2505 


2.77 


2.141 


20.98 ± 0.05 


-1.36 + 0.06 


Zn 


105 


Sinil808 


Q 0528-250 


LBQS 0528-2505 


2.77 


2.811 


21.35 + 0.07 


-0.91 ± 0.07 


Zn 


304 


Snil253 


Q 055 1-366 




2.32 


1.962 


20.70 ± 0.08 


-0.35 + 0.08 


Zn 


468 


Si nil 808 


Q0841 + 129 




2.50 


1.864 


21.00 + 0.10 


-1.51+0.11 


S 


32 


Snil259 


Q 0841+129 




2.50 


2.375 


21.05 + 0.10 


-1.59 + 0.10 


Zn 


37 


Fe nil 125 


Q0841 + 129 




2.50 


2.476 


20.80 + 0.10 


-1.60 + 0.10 


Zn 


30 


S nil 259 


Q 0913+072 


LBQS 0913+0715 


2.78 


2.618 


20.35 + 0.10 


-2.59 + 0.10 


Si 


22 


Sinil526 


Q 1036-229 


CTS 460 


3.13 


2.778 


20.93 ± 0.05 


-1.36 + 0.05 


S 


80 


FenilOSl 


Q 1108-077 


BRI 1108-0747 


3.92 


3.482 


19.95 + 0.07 


-1.63 + 0.07 


Si 


49 


Sinil526 


Q 1108-077 


BRIl 108-0747 


3.92 


3.608 


20.37 + 0.07 


-1.59 + 0.07 


Si 


31 


Fenil608 


Qllll-152 


CTS 298 


3.37 


3.266 


21.30 + 0.05 


-1.65 + 0.11 


Zn 


140 


Sinil020 


Q1117-134 


BR 1117-1329 


3.96 


3.350 


20.95 + 0.10 


-1.41+0.11 


Zn 


43/45 


Si nil 808/Fe 11^1608 


Q 1157+014 




1.99 


1.944 


21.80 + 0.10 


-1.44 + 0.10 


Zn 


89 


Feni2260 


Q 1209+093 


LBQS 1209+0919 


3.30 


2.584 


21.40 + 0.10 


-1.01+0.10 


Zn 


214 


Si nil 808 


Q 1210+175 


LBQS 1210+1731 


2.54 


1.892 


20.70 ± 0.08 


-0.93 + 0.08 


S 


62 


Snil253 


Q 1223+178 


LBQS 1223+1753 


2.94 


2.466 


21.40 + 0.10 


-1.63 + 0.10 


Zn 


91 


Sinil808 


Q 1232+082 


LBQS 1232+0815 


2.57 


2.338 


20.90 ± 0.08 


-1.43 + 0.08 


S 


85 


Si nil 808 


Q 1331+170 




2.08 


1.776 


21.15 + 0.07 


-1.28 + 0.08 


Zn 


75 


Feni2374 


Q 1337+113 


LBQS 1337+1121 


2.92 


2.508 


20.12 + 0.05 


-1.81+0.06 


Si 


32 


Feni2344 


Q 1337+113 


LBQS 1337+1121 


2.92 


2.796 


21.00 + 0.08 


-1.86 + 0.09 


Si 


42 


Fenil608 


Q 1340-136 


CTS 325 


3.20 


3.118 


20.05 ± 0.08 


-1.42 + 0.08 


S 


153 


Fe nil 608 


Q 1409+095 


LBQS 1409+0930 


2.85 


2.019 


20.65 + 0.10 


-1.62 + 0.16 


Zn 


39 


Feni2374 


Q 1409+095 


LBQS 1409+0930 


2.85 


2.456 


20.53 ± 0.08 


-2.06 ± 0.08 


Si 


69 


Sinil526 


Q 1441+276 


PSS J 1443+2724 


4.42 


4.224 


20.95 + 0.10 


-0.63 + 0.10 


S 


130 


Fe nil 081 


Q 1444+014 


LBQS 1444+0126 


2.21 


2.087 


20.25 + 0.07 


-0.80 + 0.09 


Zn 


294 


Fenil608 


Q 1451+123 


LBQS 1451+1223 


3.25 


2.255 


20.35 + 0.10 


-1.07 + 0.15 


Zn 


103 


Feni2374 


Q 1451+123 


LBQS 1451 + 1223 


3.25 


2.469 


20.40 + 0.10 


-2.27 + 0.13 


Si 


35 


Feni2382 


Q 1451+123 


LBQS 1451 + 1223 


3.25 


3.171 


20.20 + 0.20 


-2.10 + 0.21 


Si 


45 


Sinil526 


Q 2059-360 




3.09 


2.507 


20.29 ± 0.07 


-1.85 + 0.20 


S 


25 


Feni2344 


Q 2059-360 




3.09 


3.083 


20.98 ± 0.08 


-1.77 + 0.09 


S 


44 


Fe nil 121 


Q2116-358 


HES 21 16-3550 


2.34 


1.996 


20.10 + 0.07 


-0.34 + 0.11 


Zn 


177 


Fenil608 
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Table 1. continued 



Quasar 


Other name 


z 


em 


, 1 


logA^(Hi) 


[X/H]^ 




X 


Ay 

(km s ' ) 


Selected 
transition lines 


Q 2138-444 


LBQS 2138-4427 


3, 


.17 


2, 


.383 


20.60 + 0.05 


-1, 


.15 + 0. 


.10 


Zn 


131 


Fe 11/12374 


Q2138-444 


LBQS 2138-4427 


3, 


.17 


2, 


.852 


20.98 + 0.05 


-1, 


.74 + 0, 


.05 


Zn 


58/62 


Si n^l808/Fe 11/12374 


Q2152+137 


PSS J2155+1358 


4, 


.26 


3, 


.316 


20.50 + 0.15 


-1, 


.37 + 0. 


.15 


Si 


74 


Fen/11608 


Q 2206-199 


LBQS 2206-1958 


2, 


.56 


1, 


.921 


20.67 ± 0.05 


-0 


.54 + 0. 


.05 


Zn 


136 


Si 11/11808 


Q 2206-199 


LBQS 2206-1958 


2, 


.56 


2, 


.076 


20.44 ± 0.05 


-2, 


.32 + 0. 


.05 


Si 


20 


Si 11/11304 


Q 2230+025 


LBQS 2230+0232 


2, 


.15 


1, 


.864 


20.90 + 0.10 


-0 


.81+0, 


.10 


S 


148 


S 11/11253 


Q 223 1-002 


LBQS 2231-0015 


3, 


.02 


2, 


.066 


20.55 ± 0.07 


-0 


.65 + 0, 


.07 


S 


145 


S 11/11253 


Q 2243-605 


HES 2243-6031 


3, 


.01 


2, 


.331 


20.65 ± 0.05 


-0 


.85 + 0, 


.05 


Zn 


1 

173 


„ 11 OAO 

Si 11/11808 


Q 2332-094 


FIRST J 2334-0908 


3, 


.32 


2, 


.287 


20.07 ± 0.07 


-0 


.33 + 0. 


.07 


Zn 


43 


Si 11/11808 


Q 2332-094 


FIRST J 2334-0908 


3, 


.32 


3, 


.057 


20.50 ± 0.07 


-1 


.33 + 0. 


.08 


S 


111 


Fen/11608 


Q2343+125 




2, 


.51 


2, 


.431 


20.40 ± 0.07 


-0 


.92 + 0. 


.07 


Zn 


289 


Fen/11608 


Q 2344+125 




2, 


.76 


2, 


.538 


20.50 + 0.10 


-1 


.81+0. 


.10 


Si 


69 


Fen/11608 


Q 2348-011 


LBQS 2348-0108 


3, 


.01 


2, 


.426 


20.50 + 0.10 


-0 


.64 + 0. 


.10 


S 


248 


Fen/11608 


Q 2348-011 


LBQS 2348-0108 


3, 


.01 


2, 


.615 


21.30 + 0.08 


-2, 


.02 + 0. 


.08 


Si 


95/105 


Si ii/11808/Si 11/11 526 


Q 2348-147 




2, 


.94 


2, 


.279 


20.63 + 0.05 


-2, 


.10 + 0. 


.08 


S 


55 


Fe n/12586 


Q 2359-022 


LBQS 2359-0216 


2, 


.81 


2, 


.095 


20.65 + 0.10 


-0 


.84 + 0. 


.13 


Zn 


146 


Fen/11608 


Q 2359-022 


LBQS 2359-0216 


2, 


.81 


2, 


.154 


20.30 + 0.10 


-1 


.62 + 0, 


.10 


Si 


67 


Sin/11526 



* Total neutral hydro gen column density and metalUcity measurements are from this work, except for the systems at Zab,s = 3.390 

toward 0000-263 ("Molaro et al."200lV Zab.s = 3.025 and 2.087 toward, respectively, Q 0347-383 and Q 1444+014 

(Ledo u x et al. 2003) . Zabs = 1 -962 toward Q 055 1 -366 cLedoux et al.„200Z) . and Zabs = 4.224 toward Q 1441+276 jLedoux et alJ 
l2006l) . 

Quoted uncertainties are Icr standard deviations. 

' There are a few DLA systems with c|zem - Zabsl/(1 + ^abs) < 5000 km s"' . These systems were not considered during the 
correlation analyzes. However, their properties (metallicity, A^(H i), ionisation level) do not resemble those of the associated 
systems. They are not associated with the central engine of the quasar, nor ejected by the quasar, but rather associated with 
dense gas in its surroundings (see text). 

^ Average gaseous metalHcities relative to Solar, [X/H] = log[A^(X)/A^(H)] - log[A^(X)/A^(H)lo, calculated using Solar 
abundances listed in Morton (2003), which are based on meteoritic data from Grevesse & Sauval (.200 2). The reference element 
was taken to be X=Zn when Zn ii is detected, or else either S or Si was used. 

^ Transition lines used to determine the velocity widths of low-ionisation line profiles. Whenever none of the observed lines 
satisfies the adopted criterion (see Sect.|2}, we used the mean value of the velocity widths calculated from (i) a Une slightly 
more saturated, and (ii) a line slightly weaker, than what the criterion specifies. 



